Non-equilibrium Thermodynamics of Spacetime 



Christopher Eling 1 , Raf Guedens 2 , and Ted Jacobson 1 
1 Department of Physics, University of Maryland 
College Park, MB 20742-4111 USA and 
2 High Energy and Elementary Particle Division, University of Crete, 
P.O. Box 2208, GR-710 03 Heraklion, Crete, GREECE 



O 
O 

(N 



O . 
O ■ 
(N ■ 

o : 

O ■ 

"o ! 

i < 
5h . 

bo 



It has previously been shown that the Einstein equation can be derived from the requirement that 
the Clausius relation dS = 8Q/T hold for all local acceleration horizons through each spacetime 
point, where dS is one quarter the horizon area change in Planck units, and 5Q and T are the 
energy flux across the horizon and Unruh temperature seen by an accelerating observer just inside 
the horizon. Here we show that a curvature correction to the entropy that is polynomial in the Ricci 
scalar requires a non-equilibrium treatment. The corresponding field equation is derived from the 
entropy balance relation dS = SQ/T + diS, where cUS is a bulk viscosity entropy production term 
that we determine by imposing energy-momentum conservation. Entropy production can also be 
included in pure Einstein theory by allowing for shear viscosity of the horizon. 



I. INTRODUCTION 

A profound connection between gravitation and ther- 
modynamics was first suggested by the discovery in the 
1970's of black hole entropy^, the four laws of classi- 
cal black hole mechanics Q, and Hawking radiation 
But it is rather mysterious that the Einstein equation, a 
hyperbolic second order partial differential equation for 
the spacetime metric, has a predisposition to thermody- 
namic behavior. A decade ago one of us proposed to 
explain this connection by reversing the logic, using the 
assumed proportionality of entropy and horizon area for 
all local acceleration horizons (called there local Rindlcr 
horizons) to derive the Einstein equation as an equilib- 
rium equation of state. This derivation suggests the idea 
that gravitation on a macroscopic scale is a manifestation 
of thermodynamics of the vacuum. 

Does this thermodynamic derivation of the Einstein 
equation indicate something deep, or is it a case of "as- 
suming the answer" , or a superficial consequence of the 
assumptions, or perhaps just an accident? We address 
this question here by investigating whether and how the 
derivation can be generalized to allow for the higher 
curvature terms expected in the field equation on the 
grounds of effective field theory One might guess that 
this could be done simply by allowing for curvature to 
enter the ansatz for horizon entropy. The purpose of this 
article is to investigate whether this is so. We find that 
entropy dependence on the Ricci scalar can indeed be 
accommodated, but it requires a change of setting from 
equilibrium to non-equilibrium thermodynamics. 

We begin by reviewing the two hypotheses on which 
the derivation of the Einstein equation of state given in 
Ref. is based. Next we present that derivation, and 
then generalize it to allow for dependence of the entropy 
on the Ricci scalar. We end with several comments. 

The motivating idea is that the origin of the ther- 
modynamic behavior of black holes is to be traced to 
the thermal nature of the Minkowski vacuum. The vac- 
uum is the ground state of the generator of time trans- 



lations, but it is thermal with respect to any genera- 
tor of Lorentz boosts. More precisely, restricted to the 
"Rindler wedge" x > \t\ (using Minkowski coordinates) 
the vacuum density matrix for a relativistic quantum field 
has the form of the canonical ensemble (Gibbs state) 
p = Z^ 1 exp(— Hb/T), where Hb is the boost hamil- 
tonian and the "temperature" is T — h/2-K. 0,0- This 
"temperature" does not have dimensions of energy, be- 
cause the boost Hb generates translations of a dimen- 
sionless hyperbolic angle, rather than of a time. When re- 
scalcd to generate proper time translations on the world- 
line of a uniformly accelerated observer this becomes the 
Unruh temperature Tjj = ha/lit , where a is the ac- 
celeration. 

The past boundary x = t < of the Rindler wedge 
is a lightlike hyperplane in spacetime, forming a causal 
boundary or horizon for the past of the "bifurcation 
plane" x = t = 0. The state behind the horizon is hid- 
den to outside observers with access only to a spatial 
slice bounded by this plane. For such observers the rel- 
evant state is the density matrix p. The entropy of p is 
infinite in standard quantum field theory, but when UV 
regulated it is proportional to the area of the bifurca- 
tion plane and depends on the number and nature of the 
quantum fields @ . The similarity to Bekenstein-Hawking 
black hole entropy, which is universally given by the hori- 
zon area divided by AhG, motivates our first hypothesis: 
we suppose that in quantum spacetime — whatever that 
is — there is a universal entropy density a per unit horizon 
area. We imagine that a indeed depends on the number 
and nature of quantum fields, if any such freedom exists 
in the underlying fundamental theory. 

Our second hypothesis concerns the relation between 
this entropy and the flux of boost energy across the 
Rindler horizon. Under a small perturbation of any 
Gibbs state at temperature T, the variation of the en- 
tropy is related to the variation of the mean energy by 
SS = 5{E)/T. When the energy is "heat" this identity 
expresses the Clausius relation dS = SQ/T. The de- 
tailed nature of energy that flows into a Rindler wedge 
cannot be examined by the outside observers, hence it 
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can be considered by them as heat that has flowed into 
the thermal system behind the horizon. This motivates 
our second hypothesis: the Clausius relation holds for all 
local causal horizons (defined precisely below) in a suffi- 
ciently small neighborhood of the bifurcation plane, with 
5Q interpreted as the mean flux of boost energy across 
the horizon and dS as the difference of area for the wedge 
with and without the boost energy flux. This hypothesis 
is inconsistent with the assumption of a fixed, flat space- 
time, since a Rindler horizon has fixed area, but it is 
consistent with the horizon-focusing effects of spacetime 
curvature provided the Einstein equation holds p| . 

We define a local causal horizon at a point p as fol- 
lows: choose a spacelike 2-surface patch B including p, 
and choose one side of the boundary of the past of B. 
Near p this boundary is a congruence of null geodesies 
orthogonal to B. These comprise the horizon. We choose 
B such that this congruence has vanishing expansion 9 
and shear a ab at p, which corresponds to an equilibrium 
state at p. We consider transitions that terminate in this 
equilibrium state. 

To define the heat flux and temperature we employ an 
approximate boost Killing vector field \ a that vanishes 
at p and whose flow leaves invariant the tangent plane 
Bp to B at p. That is, the covariant derivative Xa-,b is 
a timelike antisymmetric tensor orthogonal to B p at p. 
We normalize x° by Xa;bX a ' b = — 2, like the usual boost 
Killing vector xdt + td x in Minkowski spacetime. These 
conditions at p would determine a unique boost Killing 
vector up to sign in flat spacetime. In a generic curved 
spacetime there are no Killing vectors, so the most we 
can do is solve Killing's equation Xa:b + Xb-a = with 
this "initial data" out to some order in the neighborhood 
of p. In Ricmann normal coordinates {x a } based at p 
these initial conditions determine the zeroth and first or- 
der parts of x°j while Killing's equation imposes that 
the second order part vanishes. Generally the equation 
cannot be satisfied at third order. 

We choose the direction of x a to be future pointing on 
the causal horizon. In flat spacetime, it would then be 
related on the horizon generator through p to the affinely 
parametrized horizon tangent vector k a via x a = — Xk a , 
where A is an affine parameter that is negative and in- 
creasing along the horizon and vanishes at p. Up to the 
0(x 3 ) ambiguities in x a , the same is true in the curved 
spacetime. 

As motivated above, we define the heat as the mean 
flux of the boost energy current of matter across the hori- 
zon, 
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jT abX a dY>\ 



(1) 



where T ab is the expectation value of the matter stress 
tensor. (We omit the bracket notation for brevity.) This 
and all subsequent integrals are taken over a short seg- 
ment of a thin pencil of horizon generators centered 
on the one that terminates at p. Using the relation 
X a = —Xk a (which holds on the generator through p up 



to order 0(x 3 )) and T — h/2w we thus have 

^ = (2n/h) J T ab k a k\~X)dXd 2 A. (2) 

The entropy change SS = a 5 A is determined by the 
area change of the horizon, 



5 A = / 6dXd 2 A 1 



(3) 



where 9 = d(\nd 2 A)/dX is the expansion of the congru- 
ence of null geodesies generating the horizon. Using the 
Raychaudhuri equation 
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and the assumed vanishing of 9 and a ab at p we have 

9 = -XR ab k a k b + 0(\ 2 ). (5) 
The entropy change is thus given to lowest order in A by 



SS = 



Rabk a k b (-X)dXd 2 A. 



(6) 



If we now require that SS = 5Q/T hold for all local 
Rindler horizons through all points p, we infer that the 
integrands of (J2J and © must match for all null vectors 
k a . The integrands are both first order in A, and equality 
of the coefficients of A implies the relation 



Rab + ^9ab = {2ir/Ha)T ab 



(7) 



where $ is a so far undetermined function. To deter- 
mine <E> we require that the matter stress tensor is di- 
vergence free, corresponding to the usual local conser- 
vation of matter energy. Taking the divergence of both 
sides of and using the contracted Bianchi identity 
Rab' a = ^R.b we then find that $ = — ^R — A, where 
A is a constant. Therefore {7J) corresponds to the Ein- 
stein equation R ab — ^Rg a b — J^9ab = 8nGT ab with (un- 
determined) cosmological constant A and with Newton's 
constant determined by the universal entropy density a, 



G = {Aha) 



(8) 



Note that the entropy area density is thus universally 
(4fi,G) _1 , no matter what is the nature and number 
of quantum fields, in agreement with the Bekenstein- 
Hawking black hole entropy. This completes our review 
of Ref. 0. 

We now seek the thermodynamic equation of state if 
the entropy density is taken to be a times a function 
f(R) = 1 + 0(R) of the Ricci scalar, instead of a constant 
a as before. In this case the entropy change is given in 
analogy with @ by 



(9) 



5S = a (9f + /) dXd A, 
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where the overdot signifies derivative with respect to A. 
If the expansion 9 vanishes at p then the integrand of 
at p is / = f'(R)k a R ya , which is generally non-zero. 
This cannot match the SQ/T integrand which is of order 
A. Thus the Clausius relation implies that 9(p) must be 
non- vanishing, so as to cancel off the derivative of f |T^. 
That is, we must have 

(9f + f)(p) = 0. (10) 

This means that the causal horizon must be defined as 
the boundary of the past of a 2-surface B that is warped 
at p such that (|10l) is satisfied. This does not coincide 
with the approximate Killing horizon as closely as when 
the expansion vanishes. Nevertheless, the approximate 
Killing vector is still related to the tangent of the horizon 
generator through p by \ a — — Afc a , up to the 0(x 3 ) 
ambiguity of \ a ■ 

Since the area of the horizon is changing at p, it seems 
at first that the "system" does not approach an equi- 
librium state at p, even though the entropy is instanta- 
neously stationary there. However, the relevant notion 
of time here is the Killing flow. The relation between the 
affine parameter A and Killing parameter v on a Rindler 
horizon is A = — exp(— v), so the point p occurs at infi- 
nite Killing time. Even if 9(p) ^ 0, the rate of change of 
area with respect to Killing time vanishes like ~ exp(— v) 
as p is approached. This can be considered an approach 
to equilibrium. In the special case 9(p) — considered 
previously, the expansion vanishes at twice this rate, i.e. 
as ~ exp(— 2v). The slower decay rate in the general case 
suggests that equilibrium thermodynamics may not ap- 
ply, so that the Clausius relation may not hold. Instead 
we may have dS > 5Q/T, or more precisely the entropy 
balance relation 

dS = SQ/T + d l S, (11) 

where diS is the entropy developed internally in the sys- 
tem as a result of being out of equilibrium 0. The 
entropy production rate vanishes at p, since that is an 
equilibrium point, so we expect the rate is of order A. 

To extract the O(A) term in the integrand of we 
differentiate with respect to A and use l|10fl . 

^(0/+/)U=o = 0f-r 1 / 2 +/' (12) 

Using the Raychaudhuri equation JJJ and the geodesic 
equation k a k b - a = this takes the form 

~k a k\fR ab - f, ab + r\f, a ib) - ye 2 . (13) 

Were there no entropy production diS, the Clausius re- 
lation would imply that (|13fl must be equal to the coef- 
ficient of A in the heat flux integrand of J2J) for all null 
vectors k a . It would follow that in place of J7J) we have 

fRab - fab + f/" 1 /^ f,b + ^9ab = (2ir/ha)T ab , (14) 

where (|l()|l has been used to re-express the 9 2 term and 
$ is a so far undetermined function. We now show that 
this is inconsistent with energy conservation. 



We require, as before, that the matter stress tensor is 
divergence free, so the divergence of the left hand side 
of (jl4|l must vanish. Using the contracted Bianchi iden- 
tity, the commutator of covariant derivatives 2v c .[ ab ] = 
R a bd c v d , and defining L by / = dC/dR, we find 

(fRab - Ub)^ = {\c - □/),„. (15) 

Thus we must have 

# = D/-^C-e, (16) 

where the gradient of O matches the divergence of the 
remaining term in i|14|) , 

e, 6 - (^jiai b y a . (i?) 

This reveals a contradiction, however, since the right 
hand side of l|17|) is generally not the gradient of a scalar. 

We propose that this contradiction with energy- 
momentum conservation is resolved by the entropy pro- 
duction term diS in (|llfl . Examination of l|13l) shows 
that the problematic term would be canceled if we set 

d t S = J ad\d 2 A (18) 

with entropy production density 

a=-\a]- 1 j 2 \ = -\aJ9 2 \ (19) 

(using iJTDJ). In terms of the expansion with respect to 
Killing parameter, 9 — 9(d\/dv), we have 

adX^^af9 2 dv. (20) 

This is just like the entropy production term for a fluid at 
temperature T due to a bulk viscosity rj — (3/2)a/T 0. 
Putting T equal to the boost temperature H/2tv yields 
t] = 3haf/An. With this for diS, equation {TTJ at 0(A) 
implies the equation of state 

fRab - f-ab + (□/ - \c)g ab = {2n/ha)T ab . (21) 

This coincides with the equation of motion arising from 
the Lagrangian (ha./Aii)C(R) for which the entropy den- 
sity of a stationary black hole horizon is af(R) [l(J]. The 
thermodynamic equation of state is therefore again con- 
sistent with the Lagrangian field equation, as in the pure 
GR case. We conclude with a number of remarks. 

1. Given an entropy functional of the macroscopic vari- 
ables of an ordinary thermodynamic system, one can nor- 
mally derive the equation of state from the Clausius re- 
lation together with the first law of thermodynamics. In 
our case the first law was not explicitly invoked. How- 
ever, to fix the trace part of the field equation we required 
that the energy-momentum tensor of matter is divergence 
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free, which expresses the local conservation of matter en- 
ergy and is hence a form of the first law. 

2. It is common in near-equilibrium thermodynamics 
to deduce the general form of entropy production terms 
quadratic in the gradients of state variables, but the coef- 
ficients of those terms are phenomenological and depend 
on details of the microphysics 9] . It may therefore seem 
puzzling that here we deduce also the bulk viscosity r\. 
However, r\ is precisely 3h/Air times the entropy density 
a/, and a is purely phenomenological in our derivation. 
The only puzzle is therefore the simple relation between 
entropy density and bulk viscosity. Presumably general 
covariance lies at the root of this. 

3. In pure GR the bulk viscosity appears to become 
3Sa/47r, but this is not correct since 6(p) — 0, so there 
is no O(A) entropy production term in the integrand of 

By contrast, for globally defined black hole horizons 
it has been shown that the bulk viscosity is negative 
and equal to — l/lQnG = — ha/ An. 

4. In the presence of curvature contributions to the 
entropy, we had no choice but to use a non-equilibrium 
entropy balance relation, due to the nonzero expansion 
at p I|1U|) . Could a non-equilibrium description be used 
voluntarily also for pure GR, where the entropy is just 
the area? We cannot allow for expansion at p, since with- 
out the derivative of / to balance the expansion, there 
would be a first order term in the entropy change 5S not 
matched by the heat flux term. But how about allow- 
ing for shear at pi We initially set the shear to zero 
on the grounds that it was required by equilibrium at p. 
This is not valid however, since the shear defined with 
respect to Killing parameter would in any case vanish. 
The consequence of nonzero shear at p is to introduce 
into JSJ a term of the form — \a a bO' ah ■ This additional 
term can be written in terms of derivatives of k a , which 
can be independently chosen at p. The k a parts of the 
Clausius relation imply the Einstein equation as above, 
but the dk a part is satisfied only if the shear vanishes. 
However, nonzero shear at p is allowed if we include an 
internal entropy production term aa a bcr ab correspond- 
ing at temperature T(— hj2ii) to a shear viscosity 
Ta/2 = ha/4ir = l/16irG, just as for a black hole hori- 
zon [m. 



5. Under what conditions are the higher curvature 
terms meaningful in the thermodynamic interpretation? 
Let the entropy expansion coefficients /3j be defined by 
f(R) = 1 + (3iR + /3 2 R 2 + ■■■ , and let the curvature 
length scale at p be set by L c . Then f, a f,b in lfl4l is 
0(/3i/L 2 ) smaller than / )0 & and the 0(R 2 ) part of fR a b, 
which are smaller than the Einstein term by the same 
factor. Our assumption that the entropy balance rela- 
tion holds over an arbitrary small patch of local causal 
horizon is questionable due to quantum fluctuations at 
the Planck scale. If we assume only that the relation 
holds over a region of some minimum coordinate size e 
in Riemann normal coordinates, the definition of heat 
flux will be fuzzy. How fuzzy? The Killing vector is am- 
biguous at 0(x 3 ), while the lowest order part is O(x). 
The corresponding fuzziness in the heat flux is therefore 
of relative size 0{e 2 /L 2 ). The ratio of this to the rela- 
tive suppression of the L~ and L~ 6 terms in the field 
equation is e 2 //3i and (eL c //3i) 2 respectively. If f3i is suf- 
ficiently large (/3j 3> eL c ) all the curvature corrections 
will be meaningful as long as e < L c . This is fine for a 
"theoretical laboratory," however, dimensional analysis 
suggests that in nature we have (3\ ~ e 2 ~ -^pianck- ^ 
this is the case, then the fuzz is always comparable to 
the L~ 4 terms and much larger then the L~ 6 ones. This 
conclusion could be evaded if it were somehow correct to 
define the heat flux using the well-defined horizon gener- 
ating vector —Xk a rather than the fuzzy Killing vector. 

6. Having looked at the simplest case it would be inter- 
esting to allow for other types of curvature contributions 
that can arise. One might imagine that results similar 
to those found here would apply. However, since the 
entropy can involve not just the curvature scalar but dif- 
ferent projections of curvature into the horizon or normal 
to it, this is certainly not clear. 
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